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Abstract —The two order partial differential equation in the 
image edge detection, has good edge holding ability, but there 
will be a ’’piecewise constant" effect, leading to false edges 
produced. Partial differential equation of higher order is very 
good to overcome the "piecewise constant" effect, but the ability 
of keeping edge drop, resulting in blurred edge. Therefore, this 
paper proposes a model for edge detection fusion of two order 
partial differential equation with four order partial differential 
equations, gives the Euler Lagrange equation and the gradient 
descent flow, and introduces the numerical discretization 
method.Finally a numerical validation is performed by the edge 
detection of objective evaluation index. Experiments show 
that, this method is very effectively in the image edge 
extraction, and overcome "piecewise constant" effect, edge 
blur and other shortcomings, improving the extraction results. 

Index Terms —the two order PDE, four order PDE, edge 
detection, evaluation index 


I. Introduction 

Image edge detection, which is a key technology in the 
machine vision field [1], can make necessary early-stage 
preparation for later image processing such as image 
segmentation [2], image inpainting [3], and pattern 
recognition [4]. The traditional second-order partial 
differential equation models include the models based on heat 
diffusion equation (such as heat equation diffusion model [5], 
P-M diffusion model [6] and Catte diffusion model [7]) and 
the second-order models based on total variation (such as the 
classical TV model [8], General TV model [9] and 
self-adaption model [10]). The second-order partial 
differential equation model is able to preserve the edges well 
during the edge extraction process. However, it will generate 
a “piecewise constant effect” during the image smoothing 
process, result in the existence of false edge and, thus, affect 
the extraction effect. The four-order differential equation 
models arisen in recent years include the Y-K model [11] and 
the improved Y-K model [12]. They can successfully 
overcome the shortcomings caused by second-order partial 
differential models but will cause the shortcomings such as 
“speckle effect” [13] edge blur [14] in images and inaccurate 
edge extraction and positioning during the smoothing process. 

After the analysis on the advantages and disadvantages 
of the traditional second-order and four-order partial 
differential equations during the edge extraction process, an 
improved high-order partial differential model was obtained 
through the integration between the second-order and 
four-order models. Its gradient descent flow was given 
through the Euler-Langrange equation. Finally, the 
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experimental verification was conducted through the 
numerical discretization method and in combination with the 
objective evaluation indexes of edge detection. The result 
shows that the model is improved to some degree when it is 
compared with the traditional second-order and four-order 
partial differential models based on the objective evaluation 
indexes of edge detection. 


II. New model of improved high-order partial 

DIFFERENTIAL EQUATION 
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image smoothness. In this paper, the new model integrating 
the second-order partial differential equation and the 
four-order partial differential equation is used namely, 
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The operator is used to describe the image smoothness. 
Then, the energy functional obtained by the new model in the 
condition of noise variance constraint is: 
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In the formula above, the second item, which is a fidelity 
term, is used to approach the similarity between the processed 

image and the original one. is the regulation coefficient. 
Solving the Euler-Lagrange equation of the energy functional 
is to solve the minimum value of the functional. Change 
formula (2) into: 
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In the formula above, the integrand is: 
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Then, the first-order partial derivative and the 
second-order partial derivative in formula (4) is: 
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The Euler-Lagrange equation in formula (3) is: 
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Then, formula (7) can be simplified into: 
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Continue to simplify formula (9): 
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Therefore, the gradient descent flow of formula (10) is: 
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In order to weaken the influence of noise, Gaussian filter 
can be conducted before calculation of D, namely, 
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Then, formula (11) is changed into: 
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The numerical solution of image smoothing can be 
obtained through numerical discretization of formula (13). 

III. Numerical discretization method 

In this paper, the following difference scheme is used to 
discretize formula (13). 
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The corresponding first-order and second-order partial 
derivatives can be solved inside each bracket of formula (13) 
but the solving process must be very complicated. For the 
convenience of operation, we can directly use the first-order 
difference formula and second-order difference formula in 
formula (14) and formula (15) to solve the first-order and 
second-order partial derivatives outside the brace of formula 
(13). 
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Then, 
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In order to prevent a in formula (16) from equaling to 
0 to result in the failure of calculation, add a tiny amended 
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converted into 08 , namely 
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As the high-order partial derivative existing in the 
formula, which is used to describe the image smoothing, will 
certainly bring in the interference of “speckle effect”, we 
provide a discretization method [15] to get rid of the “speckle 
effect”; namely, calculate the four-directional mean m and 
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IV. Detection result and data analysis 

Test the teaching Lena test image with the 512x512 
international standard and the teaching building test image for 
crosswise comparison; then, compare the evaluation index of 
edge detection image obtained with the evaluation indexes of 
edge detection images respectively obtained from (1) heat 
equation diffusion model, (2) P-M diffusion model, (3) Catte 
diffusion model, (4) classical TV model, (5) General TV 
model and (6) self-adaption TV model among the traditional 
second-order models and the (1) classical Y-K model and (2) 
Y-K model improved based on the diffusion coefficient 
among the traditional four-order models. 

As for the hardware testing environment, computers with 
the Core i5 processor (2.4GHz) and 4GB of physical memory 
in the laboratory are chosen uniformly; as for the software 
testing environment, the Matlab2012b version is chosen 
uniformly. Let the number of procedure iterations be 200 
times, the time step be 0.01, the regularization parameter be 
0.01 and the gradient threshold be 15. 

A. Experiment I 

We used Lena as the test image. Namely, we used the 
original pollution-free Lena test image to verify the validity of 
the new model in the condition of not adding Gaussian white 
noise. The test results for the 9 types of models are shown in 
Fig. 1 below: 



(5) General TV model (6) Self-adaption TV model 



(9) Model stated in this paper 


Fig. 1 Comparison of Lena Image Edge Detection Results of 

9 Types of Models 

As is shown in the comparison of detection results, in the 
condition without noise interference, the edge detection 
effects obtained through the heat equation diffusion model 
and General TV(P=1.5) model are most unsatisfactory in the 
subjective evaluation method system; there are only a few of 
discontinuous edge pixel points and there is a serious lack of 
edge pixel points in the hair texture part of Lena. The P-M 
model and the Catte model are basically the same. The P-M 
method can be theoretically proven to be “morbid”; i.e. any 
subtle changes in the input will lead to complete changes in 
output. However, as this experiment is free of noise 
interference and the Lena images inputted in the two groups 
are totally the same, the edge extraction effects obtained from 
the detection of the two models should be the same. The 
classical TV model and the self-adaption model are generally 
the same and one cannot tell which one is better through 
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observation with naked eyes; as for the Y-K model and the 
improved Y-K model, the former one is obviously superior to 
the latter one in the number of edge pixel points, hair texture 
part and edge continuity etc. As for the new model put 
forward in this paper, even though there is an edge loss in on 
the left side of image, the capturing of weak edge, 
preservation of hair texture part and elimination of false edge 
are improved. 

In order to verify the conclusion made through the 
above-mentioned subjective evaluation method, the objective 

evaluation index Eldx 0 f numerical edge detection was 
used for verification. Besides, the evaluation index includes 
three benchmark evaluation parameters, which are 
respectively the edge detection image reconstruction 

similarity SSIM ? edge reliability Bldx anc [ edge 

continuity Cldx The evaluation results are shown in Table 
1 and Fig. 2. 

It can be seen that the comparison of numerical edge 
detection evaluation indexes verified the conclusion made 
through the above-mentioned subjective evaluation method: 

(1) The edge detection evaluation indexes of heat equation 
diffusion model and General TV (P=1.5) model are the 
smallest. 

(2) P-M model and Catte model are completely the same in 
the edge detection image reconstruction similarity index, edge 
reliability index and edge continuity index; i.e. they are totally 
the same. 

(3) As for the classical TV model and self-adaption TV 
model, if one cannot tell which detection result is better 
through observation with naked eyes, the evaluation index can 
correctly show that the classical TV model is superior to the 
self-adaption TV model. 

(4) Through the comparison of edge detection evaluation 
indexes of the classical Y-K model and the improved Y-K 
model, the latter one is far superior to the former one. 

(5) The edge detection evaluation index Eldx 0 f th e new 
model put forward in this paper is the biggest among all 
models. 

Table 1 Edge Detection Evaluation Indexes of Lena 


Image of 9 Models 


Image 

Image with noise 

Evaluation index 

SSIM 

Bldx 

Cldx 

Eldx 


Heat equation 
diffusion 

0.8245 

0.1917 

0.9171 

0.5701 


P-M 

0.8998 

0.1597 

0.9723 

0.7187 


Catte 

0.8998 

0.1597 

0.9723 

0.7187 


Classical TV 

0.8345 

0.2165 

0.9718 

0.7247 

Methods 

General TV 
(P=1.5) 

0.5795 

0.2806 

0.8426 

0.5704 

Self-adaption 
TV 

0.8348 

0.2163 

0.9378 

0.6992 


Classical Y-K 

0.8163 

0.1777 

0.9485 

0.6801 


Improved Y-K 
model 

0.9035 

0.1552 

0.9852 

0.7243 


Model stated in 
this paper 

0.8491 

0.1657 

0.9558 

0.733 


* SSIM _ e dg e detection image reconstruction similarity; 
Bldx _ edge reliability; Cldx _ e dge continuity; objective 
evaluation index of edge detection Eldx 
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Fig. 2 Curve Graphs for Edge Detection Evaluation Indexes 

of Lena Image of 9 Models 

B. Experiment II 

For the purpose of crosswise comparison, we used a 
teaching building as the test image (as shown in Fig. 3) 
without Gaussian white noise; i.e. we used the original 
pollution-free main teaching building test image to verify the 
validity of new model. The test results of the 9 models are 
shown in Fig. 4. 

It is found through the comparison of detection results 
that the edge detection effects obtained from the heat equation 
diffusion model and the General TV (P=1.5) model are most 
unsatisfactory and the edge loss is most severe in the 
subjective evaluation method system under the condition 
without noise interference. The P-M model and the Catte 
model are basically the same and one cannot tell whose effect 
is better with naked eyes. The classical TV model and the 
self-adaption model are generally the same; as for the 
classical Y-K model and the improved Y-K model, the former 
one is obviously superior to the latter one in the number of 
edge pixel points and texture preservation function. The new 
model put forward in this paper is superior to the other eight 
ones in the number of edge pixel points, texture edge 
preservation function and edge continuity. 

In order to verify the result obtained through the 
above-mentioned subjective evaluation method, we continued 
to use the objective evaluation index of numerical edge 
detection for verification, as shown in Table 2 and Fig. 5: 

It can be found that the comparison of numerical edge 
detection image evaluation indexes mentioned above can still 
lead to the conclusion made through the above-mentioned 
subjective evaluation method. 

(1) The edge detection evaluation indexes of heat 
equation diffusion model and General TV(P=1.5) are still the 
worst among all models. 

(2) The P-M model and Catte model are completely the 
same in the edge detection image reconstruction similarity 
index, edge reliability index and edge continuity index; 
namely, they are totally the same. This re-verifies that the P-M 
model and Catte model have the same edge extraction effect 
without noise interference. 

(3) As for the classical TV model and self-adaption TV 
model, if one cannot tell whose detection effect is better with 
naked eyes, the evaluation index can correctly show that the 
self-adaption model is slightly superior to the classical TV 
model. 
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(4) Through the comparison of edge detection 
evaluation indexes of classical Y-K model and improved Y-K 
model, the latter one is slightly inferior to the former one. 

(5) The new model put forward in this paper has good 
performance in the edge detection image reconstruction 

similarity index SSIM ? ec [g e reliability index Bldx anc [ 

edge continuity index Cldx j ts ec [g e detection image 

evaluation index Eldx j s i% to 4% higher than that of all 
other models. 



Fig. 3 Teaching Building Test Chart for Crosswise 

Comparison 



(1) Heat diffusion equation model (2) P-M model 



(5) General TV model (6) Self-adaption TV model 



(9) Model stated in this paper 
Fig. 4 Comparison of Teaching Building Image Edge 
Detection Results of 9 Models 


Table 2 Evaluation Indexes for Teaching Building Image 
_ Edge Detection of 9 Models _ 


Image 

Image with noise 

Evaluation index 

SSIM 

Bldx 

Cldx 

Eldx 

Methods 

Heat equation 
diffusion 

0.8656 

0.2701 

0.9755 

0.691 

P-M 

0.8837 

0.2602 

0.9951 

0.7238 

Catte 

0.7238 

0.2602 

0.9951 

0.7238 

Classical TV 

0.8746 

0.314 

0.9933 

0.7437 

General TV (P= 1.5) 

0.6698 

0.3929 

0.9622 

0.6739 

Self-adaption TV 

0.8747 

0.3149 

0.9932 

0.7439 

Classical Y-K 

0.8826 

0.2452 

0.9952 

0.7158 

Improved Y-K 

0.8724 

0.2266 

0.9931 

0.7149 

Model in this paper 

0.8852 

0.2518 

0.9949 

0.7544 


* SSIM _ e dge detection image reconstruction similarity; 
Bldx _ edge reliability; Cldx _ e dge continuity; objective 
evaluation index of edge detection Eldx 
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Fig. 5 Curve Graph for Evaluation Indexes of Building 
Image Edge Detection of 9 Models 

V. Conclusion 

After the application of the traditional second-order 
partial differential equation and the four-order partial 
differential equation in image edge detection was analyzed, a 
conclusion was made on the advantages and disadvantages of 
the second-order and high-order partial differential equations; 
then, a new edge extraction mode based on the second-order 
and high-order partial differential equations was obtained by 
preserving their advantages and eliminating their 
disadvantages, namely, by integrating the advantages of 
second-order and high-order partial differential equations. 
The model can not only effectively overcome the common 
problem of “staircase effect” of second-order partial 
differential equation but also effectively restrain the 
advantages, such as “speckle effect” and inaccurate edge 
positioning, of four-order partial differential equation. The 
Euler-Lagrange equation for the new model was solved to 
lead out the gradient descent flow, and the numerical 
discretization method was given for the purpose of the 
follow-up numerical calculation. In order to make crosswise 
comparisons, two experiments were conducted in the 
experimental process of this paper. Experiment I was to 
conduct the edge extractions for the traditional second-order 
partial differential equation, four-order partial differential 
equation and new model respectively for the Lena image 
without noise interference. Experiment II was to conduct the 
edge extractions for the traditional second-order partial 
differential equation, four-order partial differential equation 
and new model respectively for the building image. Finally, 
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the subjective and objective evaluation systems of edge 
evaluation were used to evaluate the test result. The 
evaluation result shows that the new model has good edge 
detection ability. In the future, it is required to further study 
the existence and uniqueness of solutions of second-order and 
four-order partial differential equation models as well as the 
optimal stopping time of edge extraction etc. 
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